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Abstract. IH 

This paper presents the results of spectral observa- 
tions for the largest complete sample of very luminous 
IRAS galaxies obtained to date. The sample consists of 
those 73 objects for which log(LiR./L Q ) > 11.5 (H = 
50kms -1 Mpc -1 , q = 0.5) and mag < 15.5, and was ex- 
tracted from the 2 Jy IRAS redshift catalog. All the spec- 
tra were obtained using 2.16m telescope of Beijing Astro- 
nomical Observatory during the years 1994-1996. A total 
of 123 galaxy spectra were obtained with spectral ranges 
of 4400A to 7100A and 3500A to 8100A at resolutions of 
11. 2A and 9.3A respectively. In addition to the 73 spectra 
for sample galaxies, we also present spectra for ten non- 
sample galaxies and a further 40 for the companions of 
sample galaxies. The data presented include nuclear spec- 
trum and the parameters describing the emission lines, 
absorption lines and continua as well as DSS images and 
environmental parameters. 
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1. Introduction 

The Infrared Astronomical Satellite(IRAS) all-sky survey 
provided a large database(IRAS Point Source Catalog, 
Version 2, 1988) of infrared galaxies. Based on recently 
completed redshift surveys of IRAS galaxies, statistical 
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spectroscopic analysis of very luminous IRAS galaxy sam- 
ples can now be performed. Sanders et al. (1988) have al- 
ready studied the spectra of a small sample of ten ultralu- 
minous IRAS galaxies selected from IRAS Bright Galaxy 
Survey (BGS) of Soifer et al.(1986). Recently Kim et al. 
(1995) completed a spectroscopic survey based on a large 
sample of luminous IRAS galaxies that was extracted from 
both the BGS and infrared warmed catalogs. However, the 
infrared luminosity range of Kim et al.'s sample is from 
10 10 5 L Q to 10 12 5 L Q and it therefore can not represent 
very luminous IRAS galaxies. 

This paper presents the spectroscopic data for a 
large complete sample of very luminous IRAS galaxies 
(VLIRGs), extracted from 1.936 Jy Redshift Survey of 
Strauss et al.(1992). The DSS images are also presented, 
here. A detailed analysis of the spectra and environment 
will be presented in paper II (Wu et al. 1997). 



2. Sample selection 

Strauss et al. (1990) selected 5014 objects from the IRAS 
database(IRAS Point Source Catalog, Version 2, 1988) ac- 
cording to the criteria: 
1. f 60 > 1.936 Jy 

2- ^60 > f 12^25 

a color criterion distinguishing galaxies from ob- 
jects in the Galaxy 
3. Galactic latitude |b| > 5° 



Strauss et al.(1992) then went on to publish the redshifts 
of the 2658 objects which were galaxies (here after the 2Jy 
samples) . 

Our VLIRG sample is a subset of the 2Jy sample. Con- 
sidering the observatory site, instruments and possible ob- 
servation times, we selected galaxies according to the fol- 
lowing criteria: 
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1. S > 0° 

2. log(L m /L ) > 11.5 

(H = 50kms~ 1 Mpc- 1 ,q = 0.5) 

3. mag < 15.5 

Here the magnitudes are from the 2Jy catalog. These 
criteria guarantee S/N ratios good enough for spectro- 
scopic classification. 73 of 2Jy-catalogue VLIRGs met 
these criteria. 

On account of the low S/N at 12^im and 25/zm for 
some of sources in 2 Jy sample, we used Dennefeld et al.'s 
(1986) formula for calculating the infrared flux, between 
42.5/xm and 122.5yum: 



Position Distribution 



F m = 1.75[2.55S 60 + 1.01Si 00 ]lCT 14 Wnr 



(1) 



where, S 6 o and Sioo are the fluxes at 60^m and lOOyum 
respectively. 

The infrared luminosity is therefore: 



Lir - 4^D 2 Fi R 



(2) 



for each object,where D is the distance from the Galactic 
center. 

3. The sources lists and general properties 

Table 1 lists the 83 target galaxies, of which 73 belong to 
the complete sample. Of the other ten objects in brackets, 
four of them have slightly lower infrared luminosities than 
10 11 ' 5 Lq and six are fainter than our magnitude limit ( 
mag > 15.5 ). This table presents the infrared luminosi- 
ties in unit of L©, magnitudes and redshifts. All of these 
data were derived from the 2Jy catalog, except that in 
the case of IR09517+6954 (M82). We can not use redshift 
as a distance indicator on account of the proximity. For 
this reason, we adopted the distance value given by Tully 
(1988). These objects flagged with asterisks in column 1 
are the sources which are also included in the IRAS Bright 
Galaxy Sample (S 60 > 5.24Jy) of Kim et al.(1995). A de- 
tailed comparison between these two samples will be made 
in paper II. 

The spatial distribution of the 73 sample galaxies on 
the sky is shown in Fig.l. The dotted curve is the celestial 
equator. The scarcity of objects near galactic declination 
b ~ 0° is due to the matching Zwicky catalog (1961-1968) 
which misses galaxies at low latitude. 

Fig. 2 shows the sample distribution as a function of 
redshift. The solid boxes represent the complete sample, 
while the dotted boxes include the other ten sources. It 
is clear that the redshifts concentrate in the range 0.02 
to 0.04. The median value is 0.0324 ( corresponding to 
a recession velocity 9700kms~ 1 ) which is larger than the 
value 5900kms _1 obtained by Kim et al.(1995) for their 
sample. 




Fig. 1. The distribution of very luminous IRAS galaxies in the 
sky shown in equal area projection using Galactic coordinates. 



o.i 
Redshift 



Fig. 2. Distribution of redshifts. The solid boxes represent our 
complete samples while the dotted boxes include the other ten 
galaxies. 
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Fig. 3. Distribution of infrared luminosities. The solid and 
dotted boxes represent the same galaxies as they do in Fig.2. 
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The infrared luminosity distribution is shown in Fig. 3. 
The solid and dotted boxes correspond to the same sam- 
ples as in Fig. 2. The counts decrease rapidly as the in- 
frared luminosity increases. Most of the galaxies have in- 
frared luminosities log(LiR/L Q ) between 11.5 and 12.0. 

4. Spectroscopic observations 

The spectroscopic observations were made between 
Apr.15, 1994 and Feb.15, 1996. We obtained a total of 
123 galaxy spectra. All of the observations were made with 
2.16m telescope at the Xinglong Station of Beijing Astro- 
nomical Observatory using Zeiss universal spectrograph 
with a grating of 195 A/mm dispersion. 

Before Nov.5, 1994, a Tek 512x512 CCD which covered 
a 2700A range from 4400A to 7100A at a resolution of 11.2 
A (2 pixels) was used. After that, a Tek 1024x1024 CCD 
which covered a 4600A range from 3500A to 8100A at a 
resolution of 9.3A (2 pixels) was employed. 

In most cases, slit width of about 3" was chosen to 
match the typical seeing disc at Xinglong Station, but 
occasionally, the seeing disc was smaller than 2" or larger 
than 5". Thus may affect the spectral classification. Slit 
position angles of 90° were generally used. For objects with 
two close nuclei ( separation < 2'), the slit was rotated 
such that two spectra could be obtained simultaneously. 

The seeing was about 3" to 4" on most of the observa- 
tion nights. In order to perform a relative flux calibration, 
KPNO standard stars were also observed on each night. 

Table 2 lists the observations by target objects. The 
standard IRAS name , together with one upper case let- 
ters representing the individual components ( See Fig. 9 for 
identifications) are listed in Column 1. Column 2 to Col- 
umn 10 give for each source: coordinates, observation date 
(Beijing Time), exposure time, airmass, slit width and po- 
sition angle respectively. The sources flagged with ticks in 
Column 1 1 were observed during good weather condition. 

5. Data reduction 

11 of the data reductions were performed using IRAF 
. The IRAF packages CCDRED, TWODSPEC and 
ONEDSPEC were used to reduce the long-slit spectral 
data. The CCD reductions included bias subtraction, flat- 
field correction and cosmic-ray removal. The dark counts 
were so low that their subtraction was not performed. We 
adopted Home's (1986) avariance weighed algorithm in 
order to extract the spectra. As for the aperture size, we 
adopted similar method to that of Kim et al.(1995) to 
minimize the aperture-related effects on the nuclear spec- 
tra. The apertures were varied according to the rcdshift of 
the objects so that they approached diameters of 2 kpc( 
H = 50kms _1 Mpc _1 ) for galaxies with z < 0.034, but 
were fixed at 4" for objects at larger redshifts as the see- 
ing was typically 3" to 4" and the pixel sizes was 1.3". 



The only two exceptions were for IR1 2323+ 1549 A and 
IR095 17+6954 (M82) which are quite nearby and a 2 kpc 
aperture was too large for them. 

An Fe-He-Ar lamp was used for the wavelength calibra- 
tions. More than 20 lines were used to establish the wave- 
length scale by fitting a first-order spline3 function. The 
accuracy of the wavelength calibration was better than 

lA. 

On most nights, more than two KPNO standard stars 
were used to perform the relative flux calibration. Atmo- 
spheric extinction was corrected using the mean extinction 
coefficients for the Xinglong station, that were measured 
out by BATC multi-color survey (Yan, 1995). 

The telluric O2 absorption bands at 6870A and 7620A 
did heaver effect on those emission lines at the similar ob- 
served wavelength. We therefore constructed an artificial 
- the spectrum of standard stars setting every wavelength 
to unity except these wavelength corresponding to the 
O2 bands. Division by the artificial spectrum could there- 
fore removed telluric absorption bands. In most cases, this 
technique worked well, especially for the band near 6870A. 
However in some cases, this correction seems not to have 
been very satisfactory and affected the measurement of 
the emission lines at 7620A. 

6. Spectra obtained 

Optical spectra are presented for all program objects in 
Fig. 4. They have all been corrected for telluric absorption, 
though some of corrections were not as good as we would 
have thought. 



.5DE-14 — 



.2SE-14 
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IRAF is provided by NOAO 



Fig. 5. Example of Gaussian fit for Hq+[NII] A6584 + A6548 
emission. Three narrow Gaussian components are used. The 
pluses are observed data. 



The measurements of emission line, absorption line 
and continuum strengths were performed within the IRAF 
environment using tasks (SPLOT and SPECFIT). For 
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Fig. 6. Example of multi-component fit for the galaxies with 
strong Balmer absorption. One Lorentz absorption and one 
narrow Gaussian emission components are used. 
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Fig. 7. Example of Gaussian fit for Seyfert Is. One broad and 
three narrow Gaussian emission components are used. 



-0.001 0.001 

Difference of observed and 2Jy redshift 



isolated emission lines such as [OIII]A5007, [OII]A3727 
and [OIJA6300, both Gaussian fitting and direct integral 
methods were used. For the blended lines such as Ha, 
[NII]A6548, A6584, and [SII]A6716, A6731 double lines, we 
employed a multi-Gaussian component method using task 
SPECFIT, to deblcnd them. There are three parameters 
for each Gaussian component: the central wavelength, the 
total flux and width of the emission line; and two pa- 
rameters for each continuum component: the flux and the 
slope. In order to speed up the convergence, some limit 
conditions were adopted. For example, we fixed the center 
wavelengths of several components relate to one another. 
In Fig. 5, for example, we used the simplex algorithm for 
the fitting and obtained the result via a chi-square min- 
imization process. As for the spectra with obvious H/3 
absorption indicating an underlying stellar population, a 
similar method was used. Because of the coexistence of H/3 
emission and absorption, one emission and one absorption 
component were used for the H/3 fitting. The wide absorp- 
tion wings due to stellar populations often cause the ab- 
sorption to be overestimated if Gaussian model is adopted. 
To solve this problem, we adopted the Lorentz model as 
shown in Fig. 6, and the results seem better. Some of our 
sample galaxies arc Scyfcrt-likc, and their spectra could 
not fitted well with only single Gaussian component for 
each of Ha , H/3. In those case, we combined one narrow 
(< lOOOkms" 1 ) and one broad (> lOOOkms -1 ) Gaussian 
component as shown in Fig. 7. 

The relative emission-line fluxes are listed in Table 3 
for all objects. The typical uncertainty in these measure- 
ment is about 10%. Colons (:) and semicolons (;) indicate 
values with relative uncertainties at about 30% and 50% 
respectively. For the line [OII]A3727, we could not obtain a 
value with an uncertainty of less than 20%, because it was 
at the blue end of the spectra, which can be affected by 
the low Q.E., lower S/N and poor flux calibration. In some 
cases, [OI]A6300 lines were heavily affected by the nearby 
emission lines [SIIIJA6312 which could increase the uncer- 
tainty. The telluric absorption bands 6870A and 7620A 
was also enlarge the uncertainty in the measurement of 
lines near them, despite the corrections performed. The 
double lines [SII]A6716, A6733, could sometimes be sepa- 
rated, but when this was not possible, only the combined 
values are given. 

The measured redshifts, observed Ha fluxes, equiva- 
lence widths of Ha emission lines, NalD absorption lines 
and H/3 absorption lines, and two continuum fluxes (at 
486lA and 6563A) are listed in Table 4. The typical un- 
certainty in the measured Ha fluxes was 15%, as the Ha 
lines often had to be deblended from an overlapping [Nil] 
emission line. 



Fig. 8. Distribution of redshift differences between our data 
and the 2Jy redshift catalog. 



Finally, we compare our measured redshifts with those 
of corresponding sources in the 2Jy redshift sample. The 
distribution of redshift differences is plotted in Fig. 8. The 
mean redshift difference is 0.000054 and scatter is 0.00038. 
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This means that our measurements agree well with these 
of the 2Jy catalogue. 

7. DSS images 

In order to be able to study the relationship between nu- 
clear activity and degree of interaction of VLIRG in a fu- 
ture paper in this series, we extracted the optical images 
for all our .objects from the CD-ROM version of Digital 
Sky Surveyu. The contour maps using IRAF tasks, we pre- 
sented in Fig. 9. 



Zou Z.L., Xia X.Y. and Deng Z.G., 1993, Publ. Beijing Astro- 
nomical Observatory, No. 23, 103 

Zwicky, F., et al. 1961-1968, Catalog of Galaxies and Clusters 
of Galaxies, (Pasadena:California Inst. Technology) 



8. Summary 

We have presented the observation data of optical nuclear 
spectra and DSS images of a sample of very luminous 
IRAS galaxies from 2Jy catalog. In the following paper 
(Wu et al., 1997, Paper II), we will give the results of spec- 
tral analysis and environmental study, in the same time, 
address the possible model and evolutionary sequence of 
very luminous IRAS galaxies. 
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Table 1 VLIRGs from IRAS 2Jy Catalog 



IRAS" 
(1) 


Log(Li R /L ) 
(2) 


mag 
(3) 




z 

(4) 


IRAS 
(!) 


Log(L 

( 


m/L ) 
2) 


mag 
(3) 




z 

(4) 


00189+3748 


11 


572 


15.30 





0364 


13136+6223* 


11 


937 


15 


10 





0311 


00267+3016 


11 


966 


14.80 





0504 


13183+3423* 


11 


863 


14 


80 





0230 


00509+1225 


11 


772 


14.00 





0604 


13299+1121 


11 


516 


14 


50 





0317 


01173+1405* 


11 


868 


14.90 





0312 


13362+4831* 


11 


706 


14 


10 





0275 


01324+2138 


11 


629 


15.30 





0472 


13373+0105* 


11 


701 


13 


80 





0225 


01484+2220* 


11 


851 


13.70 





0324 


13428+5608* 


12 


392 


15 


00 





0373 


01572+0009 


12 


665 


15.20 





1630 


13458+1540 


11 


821 


15 


00 





0570 


02071+3857 


11 


546 


13.00 





0179 


13496+0221 


11 


752 


15 


00 





0328 


02203+3158 


11 


837 


14.90 





0338 


13536+1836 


11 


611 


14 


80 





0497 


02222+2159 


11 


652 


14.90 





0338 


14151+2705 


11 


565 


15 


10 





0365 


02248+2621 


11 


519 


14.60 





0327 


14178+4927* 


11 


541 


15 


40 





0256 


02435+1253* 


11 


501 


15.20 





0216 


14547+2448* 


11 


897 


14 


60 





0339 


02512+1446* 


11 


780 


14.60 





0312 


14568+4504 


11 


501 


14 


60 





0357 


03117+4151 


11 


562 


14.00 





0235 


15107+0724* 


11 


525 


15 


50 





0131 


(05084+7936) 


12 


170 


15.80 





0543 


15163+4255* 


12 


072 


14 


90 





0402 


(05414+5840) 


11 


505 







0148 


15327+2340* 


12 


464 


14 


40 





0182 


(06538+4628) 


11 


490 


13.70 





0214 


15425+4114 


11 


515 


14 


20 





0317 


(07062+2041) 


11 


559 







0174 


15426+4116 


11 


546 


13 


90 





0319 


07063+2043 


11 


570 


12.60 





0173 


16104+5235* 


11 


687 


14 


50 





0294 


(07256+3355*) 


11 


467 


14.70 





0135 


16180+3753 


11 


592 


14 


30 





0307 


08354+2555* 


11 


781 


14.40 





0184 


16284+0411* 


11 


582 


14 


90 





0246 


08507+3520 


11 


811 


15.00 





0559 


16504+0228* 


12 


028 


14 


70 





0243 


(09047+1838) 


11 


490 


14.80 





0291 


16577+5900* 


11 


582 


14 


20 





0187 


09126+4432* 


11 


913 


14.90 





0393 


16589+0521 


11 


637 


15 


50 





0502 


09168+3308 


11 


725 


15.30 





0499 


17366+8646 


11 


544 


14 


60 





0264 


09320+6134* 


12 


220 


15.50 





0392 


17392+3845 


11 


554 


15 


00 





0410 


09333+4841* 


11 


523 


15.00 





0259 


17501+6825 


11 


829 


15 


20 





0512 


(09517+6954) 


10 


833 


9.57 





0009 


18525+5518 


11 


683 


15 


50 





0484 


10203+5235 


11 


620 


15.00 





0322 


18595+5048 


11 


501 


15 


10 





0271 


(10311+3507) 


12 


096 







0710 


19120+7320 


11 


624 


15 


10 





0250 


(10565+2448*) 


12 


245 


16.00 





0431 


20550+1656* 


12 


074 


15 


20 





0364 


11231+1456* 


11 


809 


15.40 





0341 


22388+3359* 


11 


531 


15 


00 





0214 


11254+1126 


11 


800 


14.80 





0410 


22501+2427 


11 


723 


15 


30 





0421 


11257+5850 


12 


040 


11.80 





0104 


23007+0836* 


11 


734 


13 


06 





0162 


11543+0124 


11 


716 


15.20 





0397 


23024+1916* 


11 


573 


15 


20 





0248 


(12112+0305) 


12 


531 







0724 


23135+2516* 


11 


730 


15 


00 





0273 


12120+6838 


12 


029 


15.40 





0608 


23254+0830* 


11 


568 


13 


60 





0290 


12251+4026 


11 


660 


15.00 





0371 


23488+1949* 


11 


528 


13 


39 





0142 


12265+0219 


12 


663 


13.07 





1583 


23488+2018* 


11 


609 


14 


90 





0179 


12323+1549 


11 


766 


15.20 





0461 


23532+2513 


11 


795 


15 


00 





0571 


12540+5708* 


12 


636 


14.10 





0418 


23594+3622 


11 


586 


15 


40 





0321 


12592+0436 


11 


787 


15.50 





0371 

















a The sources marked * in column (1) were also observed by Kim et al.(1995) 



